Direct Observation of Magnon Fractionalization in the Quantum Spin Ladder 
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We measure by inelastic neutron scattering the spin excitation spectra as a function of applied 
magnetic field in the quantum spin-ladder material (C5Hi2N)2CuBr4. Discrete magnon modes 
at low fields in the quantum disordered phase and at high fields in the saturated phase contrast 
sharply with a spinon continuum at intermediate fields characteristic of the Luttinger-liquid phase. 
By tuning the magnetic field, we drive the fractionalization of magnons into spinons and, in this 
deconfined regime, observe both commensurate and incommensurate continua. 

PACS numbers: 64.70.Tg, 75.10.Jm, 75.40.Gb, 78.70.Nx 



Gapped quantum antiferromagnets (AFs) offer some 
of the most exotic states of electronic matter ever to be 
observed, and as such have been the enduring focus of in- 
tense theoretical and experimental interest. In an applied 
magnetic field, the threefold degeneracy of the triplet ex- 
citations is lifted and the triplet gap A is closed at a crit- 
ical field B c = A/g/j,B- The properties of the states be- 
yond this quantum criticalpoint (QCP) depend strongly 
on system dimensionality and in one dimension (ID), 
where long-range order is forbidden, the field-induced 
phase is expected to be a spin Luttinger liquid (LL). In 
this state, spin-flip excitations fractionalize into spinons, 
elementary 5=1/2 entities, whose excitation spectrum 
is dramatically different from that of both ordered and 
truly quantum disordered (QD) magnets. 

Here we present the results of inelastic neutron scatter- 
ing (INS) measurements of the magnetic excitation spec- 
trum of a quantum spin ladder. By tuning the magnetic 
field, we access the full ladder phase diagram, which in- 
cludes QD, LL, and fully saturated (FM) phases. INS 
is the only experimental technique which can measure 
the full momentum- and energy-dependence of the dy- 
namical susceptibility, information which is essential to 
identify unambiguously the nature of the quasiparticles 
in these different phases, especially in the LL regime. We 
will show that both the QD phase (ladder magnetization 
to = 0) and the FM phase (to = 1) have well-defined 
magnon modes. In sharp contrast, the spectrum in the 
gapless LL state (0 < m < 1) is a continuum, whose spec- 
tral weight and incommensurate wave vector we control 
systematically. 

Previous attempts to observe the essential physics of 
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FIG. 1: Spin dynamics in the QD phase (B < B c ) at T = 50 
mK. (a-c) Triplet dispersion at B = from constant-Q scans 
along the ladder axis and in the two perpendicular directions, 
(d-e) Triplet excitations at Q — [1.5 0] in fields of and 4 
T. (f) Zeeman splitting of the triplet modes at Q — [1.5 0]. 
Solid lines are fits explained in the text. 



the field-induced LL phase in quasi-lD magnets such 
as the candidate Haldane material NDMAP [2| , effective 
S = 1 chain IPA-CuCLj 0, and possible spin-ladder 
system CuHpCl [3| have generally encountered problems 
due to additional terms in the spin Hamiltonian. These 
include single-ion anisotropy, Dzyaloshinskii-Moriya in- 
teractions, and especially interchain couplings. While a 
spinon continuum has been measured in the gapless chain 
materials Cu-benzoate @, CuPzN Q, and KCuF 3 0], 
the spin-ladder compound piperidinium copper bromide 
[(C5Hi 2 N) 2 CuBr4] [f| [§| offers the first opportunity to 
induce LL physics in a gapped system by an applied mag- 
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netic field. A number of thermodynamic measurements, 
specifically of thermal expansion , specific heat 11 1 , 
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_,, specific heat 
and magnetocaloric effect [TJ [12] , as well as by nuclear 
magnetic resonance [l3j], are consistent with predictions 
for an ideal ladder. By INS the rung and leg exchange 
parameters, J r and J;, as well as additional terms in the 
spin Hamiltonian, are determined directly. 

High-quality single crystals of (C5Di2N)2CuBr4 were 
grown from solution, and up to 10 were coaligned to ob- 
tain samples with masses of approximately 2.5 g. INS 
experiments were performed on the spectrometers IN14 
(ILL, Grenoble), FLEX (HMI, Berlin), and TASP (SINQ, 
Villigcn), using a focussing monochromator/analyzer and 
a Be filter between sample and analyzer (fixed final en- 
ergies Ef — 3.5 meV or Ef = 4.7 meV). Cryomagnets 
were used for vertical fields up to 14.8 T (_B||6~axis), and 
dilution inserts for temperatures down to 50 mK. 

INS data in the QD phase, B < B c , are summarized 
in Fig. [1] Here we focus on the ID ladder dispersion: 
sharp (resolution-limited) peaks arise from a dispersive 
triplet excitation for momentum transfers along the lad- 
der [Qh, Fig. [2(b)]. The dispersion in the perpendicular 
directions [Qk and Qi, Figs. QIa,c)] is shown to demon- 
strate the excellent one-dimensionality of the system. In 
fact we have found a systematic variation of order 30/ieV 
in the effective ladder bandwidth and subtract this 
interladder contribution to obtain the intrinsic parame- 
ters J r and J;. Figures l(d,e) show the INS intensity 
at Q = [1.5 0] (the band minimum), demonstrating a 
Zeeman splitting into three triplet components at finite 
field; Gaussian fits yield the energies shown in Fig.QJf). 

The spectrum changes dramatically above B c : we find 
a continuum of excitations extending over much of the 
Brillouin zone and up to energies of 0.8 meV. We con- 
sider first a field corresponding to half magnetic satu- 
ration (m = 0.5, Fig. [2]), where the ladder is equiva- 
lent to a gapless spin chain in zero field (below). Here 
the continuous spectrum of spinon excitations [15j is 
bounded by ei(Q h ) = fiuj(Q h ) and eJO h ) = 2?iw(%), 
where fko{Qh) = aJ\ sin (27rQ/i)| li|, 17] with a a quan- 
tum renormalization factor which is determined exactly 
from the system geometry and interaction parameters. 
In Fig. [21(a), it is clear that the commensurate m = 0.5 
spectrum is well described by such a shape. Its con- 
tinuum nature is illustrated strikingly in high-statistics 
measurements taken at B = 10.1 T: both constant-E 
[Figs.[2jb) and (d)] and constant-Q [Fig.[2](c)] scans show 
broad regions of continuous intensity, the latter extend- 
ing from 0.15 meV (lower measurement limit) to 0.8 meV. 
We stress two important points. First, these data are 
taken well inside the LL regime, at a temperature sig- 
nificantly below the LL crossover, Tll III, but above 
the boundary to 3D order induced by interladder cou- 
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continuum studied in Figs. 2 and 3 arises only from frac- 
tionalization of the lowest triplet branch of the QD phase. 
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FIG. 2: Excitation spectrum in the LL phase (Tn < T = 
250 mK < Tll) at B = 10.1 T (m « 0.5) after subtraction 
of the zero-field background, (a) Measured (left) and simu- 
lated (right) INS intensities. Solid lines mark the edges of the 
two-spinon continuum. (b,d) Constant-S scans taken along 
maxima of the transverse structure factor (insets: scan tra- 
jectories in white), (c) Constant-Q scan at Q — [0.5 0.61]. 
Black dashed lines in (b)-(d) are based on a rj-function spinon 
spectrum [solid lines in panel (a)] , red solid lines and shading 
on a full continuum calculation. 



We observe a continuous evolution of the excitation 
spectrum as the magnetic field is tuned away from m = 
0.5. Figures E3c)-(f) show constant-i? scans at B = 7.9 
T and B = 12.3 T (corresponding approximately to to — 
0.25 and to = 0.75) which are essentially identical within 
the experimental error. At E = 0.4 meV, the decrease of 
magnetic intensity compared to Fig. 2(b) between Qh = 
0.2 and 0.8 is accompanied by additional weight around 
the zone boundaries (ZBs, which we define as Qh = 
and 1). For these values of to, field-induced shifts in 
spectral weight are less pronounced at E — 0.2 meV. 
Because it is not possible by INS to follow the location of 
the zero-energy incommensurate point, instead we have 
measured the magnetic signal at the ZB as a function 
of field [Figs. [3jg,h)], finding an increase in energy and 
intensity from to = 0.5 to to = and 1. 

When the magnetic field is increased beyond a second 
QCP at B s , the spins are fully aligned. The spectrum 
becomes discrete again, as shown in Figs. HId,e), with 
elementary magnon excitations. We observe a ID dis- 
persive band [Figs. IHb,c)], whose width is very similar 
to that measured in the QD phase. The ZB excitation 
energy increases linearly with applied field [Fig. |D(a)], 
which allows the identification of B s . 

We use the magnon dispersion relations to deduce 
the exchange parameters of the system. The excita- 
tions of the spin ladder in zero field have been the 
focus of much theoretical investigation 18]. Because 
(C5Hi 2 N) 2 CuBr4 is rather "strongly coupled" (J r /Ji ~ 
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FIG. 3: Incommensurate excitations in the LL phase. (a,b) Spinon continua calculated for an S = 1/2 XXZ chain with 
anisotropy 8 — 0.5 (see text) at m = 0.25, 0.75 (a) and m = 0.5 [(b), data shown also in Fig. 2(a) convolved with instrumental 
resolution]; dashed white lines represent INS scans at E = 0.2 and 0.4 meV. (c)-(f) Constant-E scans (trajectories as in 
Fig. 2(b) inset) for B — 7.9 T and B = 12.3 T. (g) Constant-Q scans at the ZB for two chosen magnetizations m between 
and 0.5. (h) Summary of data for ZB excitation energy as a function of m. Red lines are predictions based on the XXZ chain: 
solid from panel (a) and dashed (shown for comparison) from panel (b). 



4), high-order perturbative expansions are very effective, 
and here we follow the (3D) treatment of Ref. [l9]. The 
measured triplet dispersion is dominated by the ladder 
terms: we obtain J r = 12.8(1) K and Ji = 3.2(1) K 
[black line in Fig. QJb)]. From the fits to the linear Zee- 
man splitting [solid lines in Fig. [TJf )] , the (/-factor for 
this orientation is g = 2.17(3), while B c = 6.8(1) T. The 
two approaches agree perfectly, and are consistent with 
values of g and B c determined by other techniques. 

Theoretically, the excitations of a field-polarized lad- 
der are gapped spin waves with dispersion relation 
e(Q h ) = gfJ-B(B-B s ) + Ji(l + cos(2TrQ h )) A mean- 
field treatment is exact here because all quantum fluctu- 
ations are quenched. The red lines in Figs. IUb,c) are 
fitted using this expression: when the small 3D coupling 
term is removed, we obtain J r = 13.1(1) K, J; = 3.3(1) 
K, and B s = 13.6(1) T. The fit to the ZB energy in 
Fig. [fa) yields B s = 13.6(2) T. 

The LL continuum arises from the fractionalization of 
spin-flip excitations (AS Z — 1) into two elementary and 
deconfined 5 = 1/2 objects. Indeed, in a unified descrip- 
tion of the full phase diagram, the magnon excitations in 
the QD and FM phases are bound states of these spinons, 
and the QCP at B c may be regarded as a field-driven 
spinon binding-unbinding or connnement-deconfinement 
transition. The theoretical description for a ladder in 
the strong-coupling limit may be obtained by a map- 
ping to the S = 1/2 XXZ chain with anisotropy 8=1/2 
21}. The low energy sector of the ladder [i.e. the lowest 
Zeeman-split branch in Fig. 1(f)] is governed by 



Hxxz = ^Ji(S*S? +1 +S?SV +1 +5SlS? +1 )-b eS S*, (1) 



where the effective field 6 e ff = 2b s B (J^+J^)/ 2 i s such 
that — b s < b c s < b s , with b s = | J; the saturation field 
for the XXZ chain. This effective model has spinon ex- 
citations in whose dispersion J = Ji and a = 1.299 [13], 
whence the bounds ei(Qh) and e u (Qh) shown as solid 
lines in Fig.^a). Because Si in Eq. |T]) is a composite of 
the two physical spins on each ladder rung, the total INS 
cross-section contains a rung structure factor modulating 
the contributions from longitudinal and transverse spin 
correlations 14| . All data in the LL phase were measured 
on the maxima of the transverse structure factor (insets 
in Fig. [2]) , where the longitudinal contributions are zero. 

We have calculated the transverse spin correlation 
function for all values of m following Ref. [23], and in 
Figs. 02a, b) present the results for m = 0.25, 0.75, and 
0.5. In the XXZ-chain model, the spectral intensities are 
symmetric in m about m = 0.5. Data taken at B = 10.1 
T correspond to m = 0.48(2), while fields B = 7.9 T and 
B = 12.3 T correspond respectively to m = 0.24(2) and 
to = 0.72(2). The theoretical intensities were convolved 
with the 4D instrumental resolution to obtain the global 
fit shown as red lines in Figs.[2Jb)-(d) and[3Jc)-(f), while 
the red lines in Figs. EHg,h) are obtained directly. The 
shaded bands indicate the error bar in the experimen- 
tal determination of a single constant of proportional- 
ity valid for all fields, energies, and wave vectors. Their 
width combines the statistics of all our scans with uncer- 
tainties in the exact magnetization values at the chosen 
fields and in the convolution procedure. The agreement is 
quantitatively excellent. We note in particular that only 
the smallest asymmetries between m < 0.5 and m > 0.5 
may be discerned in the data [cf. Figs. 3(c,d) and 3(e,f)], 
and that the energy-dependence of the intensity is de- 
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FIG. 4: Spin dynamics in the FM phase, (a) B-dependence of 
the ZB excitation energy. (b,c) Magnon dispersions along Qi 
and Q h . (d,e) Constant-S scans at B = 14.8 T and T = 50 
mK with Gaussian fits. Solid lines in panels (a-c) are fits to 
the dispersion, while dashed lines allow extraction of B a . 



scribed exactly [Fig. 2(c)]. 

Measurements in the QD and FM phases provide two 
independent and complementary determinations of J r 
and Ji. While the values of B c and B s deduced from 
these are marginally smaller than from thermodynamic 
studies, the latter were performed mostly on undeuter- 
ated samples and in different orientations. By exploiting 
the Q-specificity of INS, we have accounted for a 3% ef- 
fect from interladder coupling We then find that 
the intrinsic ladder leg parameters in the two regimes 
are identical within their errors, but the rung parame- 
ters are not: this apparent magnetostriction effect is of 
order 1-2%. That such a phenomenon may occur is not 
surprising in a structurally "soft" material of this na- 
ture [l(| ■ In fact this discrepancy is the sum of all addi- 
tional contributions, including any other magnetoelastic 
terms or complex spin interactions. Our magnon dis- 
persion analysis therefore quantifies the statement that 
(C5Hi 2 N) 2 CuBr4 is an excellent spin-ladder system. 

A key property of the LL spectrum is the presence 
of a zero-energy mode at an incommensurate wave vec- 
tor < Qmin < 1/2 which changes systematically with 
field. However, the spectral weight at Q m i n vanishes as 
E —> [Fig. 3(a)], precluding a direct measurement of 
the incommensurability. Instead we have presented in- 
direct confirmation of the theoretical prediction in the 
form of the field-tuned finite-E 1 spectra and the ZB en- 
ergy [Figs. 2 and 3]. The question of the evolution of 
spectral weight in the LL is of particular interest in the 
context of the commensurate 3D order ed p hase which 
emerges at sufficiently low temperature [12J, [l3| . At in- 
termediate energies, we return to the question of the sym- 
metry of the measured spectra about m = 0.5: physical 
effects arising due to departures from strict strong cou- 
pling, and from the higher triplet branches, are expected 
to cause some asymmetry in intensities, but these are 



clearly extremely small (for energies 0.5 J; < E < 2. 5 J;) 
in (C5Hi 2 N) 2 CuBr4. At high energies, it remains to ad- 
dress, both experimentally and theoretically, the nature 
of the higher spinon continua expected from the upper 
two triplet branches. 

In summary, we have performed a comprehensive INS 
investigation of the magnetic excitation spectrum in 
(C5Hi 2 N) 2 CuBr 4 , a spin-ladder compound whose energy 
scales are perfectly suited to systematic studies in labo- 
ratory fields. We observe the presence of a broad contin- 
uum of spinon excitations in the intermediate, Luttinger- 
liquid (LL) phase, which is starkly different from the 
discrete (magnon) excitations measured below the criti- 
cal field (QD phase) and above saturation (FM phase). 
From the QD and FM results, we extract the ladder 
parameters with unprecedented accuracy, demonstrating 
directly that even the sum of all other effects beyond 
the "nearly ideal" Hamiltonian of weakly-coupled lad- 
ders falls below the 2% level. The spinon continuum 
proves the occurrence of field-induced fractionalization 
into elementary 5=1/2 entities as the system enters the 
LL regime. An excellent, fully quantitative description 
of the incommensurate continua measured at all fields is 
obtained from an effective chain model for a ladder with 
the coupling ratio J r /Ji ~ 4 of (C5Hi 2 N) 2 CuBr 4 . 
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[1 
[2 
[3 
K 

[6 

ir 

[8 

[9. 
[10 

[11 

[12 
[13 
[14 
[15 
[16 
[17 

[is; 

[19 

[20 
[21 
[22 



T. 
A. 
V. 
M 
D. 
M 
B. 
B. 
B. 
T. 
Ch 



Giamarchi et. al., Nature Physics 4, 198 (2008). 
Zheludev et. al, Phys. Rev. Lett. 88, 077206 (2002). 
O. Garlea et. al, Phys. Rev. Lett. 98, 167202 (2007). 
Clemancey et al, Phys. Rev. Lett. 97, 167204 (2006). 
C. Dender et. al, Phys. Rev. Lett. 79, 1750 (1997). 

B. Stone et. al, Phys. Rev. Lett. 91, 037205 (2003). 
Lake et. al, Nature Materials 4, 329 (2005). 

R. Patyal et. al, Phys. Rev. B 41, 1657 (1990). 

C. Watson et al, Phys. Rev. Lett. 86, 5168 (2001). 
Lorenz et al, Phys. Rev. Lett. 100, 067208 (2008). 

Riiegg et al, Phys. Rev. Lett. 101, 247202 (2008). 



B. Thielemann et al. larXiv: 0809.0440v2. 

M. Klanjsek et al, Phys. Rev. Lett. 101, 137207 (2008). 

B. Thielemann et al, unpublished. 

L. D. Faddeev et al, Phys. Lett. 85A, 375 (1981). 

J. des Cloizeaux et al, Phys. Rev. 128, 2131 (1962); 

G. Muller et al, J. Phys. C: Solid State Phys. 14, 3399 

(1981). 

T. Barnes et al, Phys. Rev. B 47, 3196 (1993). 

M. Muller et al, J. Phys. Condens. Matter 12, 7633 

(2000). 

B. Normand, Acta Physica Polonica B 31, 3005 (2000). 

F. Mila, Eur. Phys. J. B 6, 201 (1998). 

J.-S. Caux et al, J. Stat. Mech. P09003 (2005). 



